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Abstract. The nearby, 30 to 40 Myr old, starburst re- 
gion of the Gould Belt has formed numerous massive stars. 
Within its ~300 pc radius, it produces core-collapse su- 
pernovae at an enhanced rate which is shown to be 75 
to 95 Myr -1 kpc -2 , i.e. 3 to 5 times higher than the lo- 
cal Galactic rate, over the past and future few million 
years. A population of persistent, but unidentified, 7-ray 
sources has been recently singled out at medium latitudes 
above 100 MeV. Their distribution across the sky is shown 
to be quite significantly and better correlated with the 
tilted Gould Belt than with other Galactic structures. As 
many as 40 ± 5 sources are statistically associated with 
the Belt at |6| > 5°. It is therefore proposed that these 
sources are part of the Belt and are relics of the Belt su- 
pernovae in the form of million-year old pulsars. Their 
presence stresses how active the local medium is, heated, 
enriched and shaped by multiple recent explosions. 

Key words: Gamma rays: observations - Stars: pulsars: 
general - Stars: supernovae: general 



1. Introduction 

Above 100 MeV, EGRET has discovered ~70 persistent, 
unidentified, sources at medium latitudes that clearly 
differ from the ~40 sources seen close to the Galactic 



plane ( Hartman ct al. 1999 ): they are significantly softer, 
fainter, and h ave a steeper logN -logS function than at 
low latitudes ( Gehrels et al. 2000 ). They have been ten- 
tatively associated with the local interstellar medium 
and the Gould Belt QGrcnicr 19911; |Gehrels et al. 200C; 



Grenier 200C), but their nature remains a mystery. The 



likely candidates present in the Galactic disc, such as pul- 
sars, supernova (SN) remnants, and OB associations, can- 
not account for so many sources off the plane. Interest- 
ingly, the starburst Gould Belt disc is tilted at ~ 20° to 
the Milky Way and close enough for weak sources to be 
detected. Proposing an origin of the sources in the Belt 
however requires one a) to confront their spatial distribu- 



tion with that of the Belt and other likely Galactic struc- 
tures given strong observational biases; b) to evaluate the 
number of likely 7-ray emitters, i.e. compact stars, formed 
in the Belt in the recent past. Both aspects are analyzed 
below before discussing the possible nature of the sources. 

2. Supernova rate in the Gould Belt 

Our Galaxy produces on average 2.5l~Q4 SNe per cen- 
tury (Fammann et al. 1994), ^85% of which arise from 
the core collapse of a massive star. This value and the 
stellar distribution in the Galaxy imply a local rate of 20 
events Myr -1 kpc -2 , in reasonable agreement with the 
29 progenitors Myr -1 kpc -2 found with masses > 8 M P) 
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within 1 kpc from the Sun ( Fammann et al. 1994 ). This 
region, however, includes the starburst Gould Belt with 
its 300 pc radius. Today, the Belt hosts 432 ±15 pro- 
genitors with masses > 8 M Q ( |Comcron et al. 1994 ) and 
their maximum lifetime will imply a crude minimum rate 
> 40 collapses Myr -1 kpc -2 in a few tens of Myr. Nu- 
merous explosions have lately occurred from the first gen- 
erations of massive stars formed in the Belt. A recent- 
past rate can be inferred from the current stellar con- 
tent of the Belt as a function of mass, for both short- 
lived and long-lived stars, given a few assumptions: 1) a 
stellar initial mass spectrum dN/dM — aA/ r_1 with a 
range of indices —2.0 < T < —1.1 as measured in nearby 
OB associations ( |Scalo 1986| ; [Massey et al. 1995| ); 2) life- 
time estimates as modelled for stars with solar metallicity 
(Z = 0.02; pchaller et al. 1992| [Meynet et al. 1994| ) and 
interpolated in mass; 3) a constant birth rate for simplic- 
ity; 4) a conservative mass threshold for collapse of 8 M Q . 
Using this formalism, star counts were integrated in var- 
ious spe ctral bands and scale d to the observed Belt star 
counts ( Comeron et al. 1994 ) to eliminate the unknown 
birth rate and amplitude a. 

Fhe significant increase in birth rate obtained for stars 
<40 Myr confirms the Belt nuclear age and the validity 
of our simple formalism. SN yields of 21.4 ± 0.8, 24.1 ± 
0.9, and 27.2 ± 1.0 per Myr were obtained for T indices 
of —2.0, —1.5, and —1.1, respectively, for a 40 Myr old 
Belt. Fhe quoted error results from the uncertainty in the 
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does influence their apparent spatial distribution. To as- 
sess their correlation with various structures, a maximum- 
likelihood test was applied that takes these biases into 



Fig. 1. visibility map, in Galactic coordinates, from the 
best fit "Iso+Bclt" model, the Belt being traced by its 
young massive stars. The colour codes the number of 
sources detectable per bin in the EGRET survey. The 67 
persistent, unidentified EGRET sources at \b\ > 2.5° are 
marked as white dots which closely follow the curved lane 
of Belt stars. 



observed star counts. These estimates decrease by 15% 
with Belt age from 50 to 30 Myr. Truncating the mass 
spectrum at 60 or 120 M Q , or choosing 10 M Q for the col- 
lapse threshold, has < 10% impact on the results. Given 
the uncertainties for the Belt age and T index (particu- 
larly at large masses), we infer a frequency of 20 to 27 
SNe per Myr and a rate of 75 to 95 Myr -1 kpc~ 2 which 
is 3 to 5 times the local Galactic one. This high rate is 
valid for the past few Myr. It is consistent with the power 
of 2.3 x 10 51 erg Myr -1 kpc~ 2 required to maintain the 



local cosmic-ray density (Blandford & Ostriker 1980) for 
a standard SN-to -cosmic-ray energy conversion efficiency 
of a few percent ( Drury et al. 1994 ). It is consistent with 
t he presence, withi n the Belt, of four 0.1-1 Myr radio loops 
( Bcrkhuijscn 1973 ) and the Local Bubble. Thirteen radio 
pulsars from the Princeton catalogue are found at high lat- 
itudes with distances <1 kpc and age <2 Myr, but they 
are too few to show a correlation with the inclined Belt 
or with the Galactic plane. The narrow radio beams from 
many more may miss the Earth. 



3. Gamma-ray sources in the Gould Belt 

The distribution on the sky of a subset of the unidenti- 
fied EGRET sources is strongly reminiscent of the tilted 
Gould Belt. The subset includes all persistent sources at 
|6| > 2.5° from the 3rd EGRET catalogue, i.e. those de- 
tected with a significance > 4er (> 5a at |6| < 10°) using 
the cumulative data from April 1991 to October 1995. 
Sources listed as likely artifacts were discarded. The 67 
sources thus selected are displayed in Fig. 1. 80% of them 
show no time variability, the others are only moderately 
variable ( Tompkins 199S ). Strong biases due to the in- 
strument exposure and the bright interstellar background 



account (Grenier 1997). An important degree of freedom 
is the unknown luminosity function dN/dL oc L Q of the 
parent sources. The flatter the function, the sharper the 
longitude and latitude profiles of any Galactic population 
because more sources remain visible to large distances. 
This effect being quite strong, a was kept as a free param- 
eter over a wide range of values: —2.5 < a < 0. EGRET 
exposure maps above 100 MeV for the 4.5-year survey 
( Hartman et al. 1999 ), the observed isotropic background 



of (2.1 ± 0.3) x 10~ 5 7 cm' 2 s" 1 sr"\ and the interstel- 



lar background model QHuntcr et al. 1997| ) of the EGRET 
survey were used to model the source visibility in 5° by 
5° bins across the sky. 

To allow for possible extragalactic sources among the 
unidentified ones, linear combinations of an isotropic and 
a Galactic component were tested (and nicknamed as in 
Table 1). The choice of Galactic distributions cover the 
various sites likely to display sources at mid latitudes: 
1) sources uniformly distributed in a spherical Galactic 
halo, 20 kpc in radius; 2) sources uniformly distributed in 
longitude in the "local Galactic disc" with any exponen- 
tial or gaussian scale height; 3) sources spread in a "thick 
Galaxy" with a gaussian radial scale length of 9.3 kpc and 
an expo nential scale heig ht of 0.4 kpc, typical of radio 
pulsars ( Lyne et al. 1985| ); 4) sources distributed in the 
interstellar mediu m as mapped by the NH=N(HI+2H2) 
column-densities ( Grenier 1997 ); 5) sources in the Gould 
Belt as traced by the column-densities of stars with spec- 
tral type < B3, fitted from the Yale Bright Star cat- 
alogue data according to the function N*(l,b) oc exp(- 
0.5*sin 2 [6 — b* (I)] /sin 2 5), with the resulting width 6 = 
12.0° ±0.4° and median latitude 6* (I) of the star distribu- 
tion displayed in Fig. 1. Distributions 4) and 5) provide 
independent ways to trace the Belt through its young stars 
as well as its clouds since the latter dominate the inter- 
stellar maps at medium latitudes. 

The reliability of the analysis was checked by applying 
it to the 67 active galactic nuclei discovered by EGRET 
and spread over the sky. For them, an isotropic distribu- 
tion indeed yields the best fit and no Galactic component, 
however large in scale height, is detected: the maximum 
log- likelihood value remains constant within 0.1 for all the 
models tested in Table 1. 67 ± 8 sources are attributed to 
the isotropic component, ± 6 to the local Gal., thick 
Gal., NH, or Belt components. The 9 ± 16 sources at- 
tributed to a halo component are not significant. The lu- 
minosity index of the parent population is a = 1.7 ± 0.3. 

Results obtained for the 67 unidentified sources are dis- 
played in Table 1. Modelled source counts in the isotropic 
and anisotropic components, -^Vj so and -/V an j respectively, 
were summed over the given latitude intervals. The prob- 
ability P Q [2|1] is determined from the log- likelihood in- 
crease between two models, with a left free. It is the 
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Fig. 2. latitude profiles of the persistent unidentified 
EGRET sources (crosses), of the best fit "Iso+Belt" model 
(thick line), and of the isotropic contribution to this 
model (dotted line). The "Iso±Local Gal." model (dashed 
line) yields a significantly poorer fit (P a [BELT|LOC] = 
3 x 10~ 6 ) with systematic deviations from the data over 
wide latitude bands. 



chance probability that a random fluctuation from model 1 
yields as good a fit as model 2. The dramatic improvement 
in the fit for any Galactic model over a pure isotropic pop- 
ulation shows that over ^50 sources have a Galactic origin. 
The fit improves very significantly from a spherical halo to 
a flatter Galactic distribution, P a [GAL | HAL] = 2x 10~ 5 , 
and even more so with the tilted geometry of the Belt 
clouds or stars, P a [BELT | HAL] = 2 x 1(T 10 . The good 
quality of the BELT fit is illustrated in Fig. 1 & 2. The 
fact that the LOC and GAL models yield equally good 
fits implies a large fraction of nearby sources: no contrast 
in longitude is detected that would require distances of 
a few kpc to the sources. Observed source counts are in- 
deed equivalent in the centre and anticentre quadrants 
in Fig. 1. Not only are the sources local, the significant 
likelihood increase between the LOC and BELT fits, with 
a chance probability P a [BELTjLOC] = 3 x 10~ 6 , gives 



strong support to their origin in the inclined Belt system. 
On the other hand, the sources at |6| > 5° were shown 
to be distinctly fainter and softer than those at lower lat- 
itudes and a subset of ^20 were pointed out along the 



Belt ( Gchrels et al. 200C ) . Together these findings provide 
compelling evidence that a distinct population of 20 to 40 
EGRET sources belong to the Gould Belt. Based on their 
spatial distribution only, they could be as numerous as 
40 ± 5 at |6| > 5°. The luminosity index of their parent 
population is a — 2.2 ± 0.3. 



4. Nature of the sources 

With luminosities L 1 of 1 to 15 x 10 26 W over 47r sr 
at 500 pc (for E~ 2 spectra), the sources cannot be un- 
resolved gas clumps irradiated by the local cosmic-ray 
flux: the required mass of - 10 4 M at 500 pc can- 
not have escaped the radio and IR surveys, even when 
considering radio beam dilution from unresolved clouds 
(Grenier 1997). Nor can they be slow (<20 km/s), old 



neutron stars, accreting gas f rom a dense c loud: they are 
at least 10 2-3 times too rare ( Grenier 1997| ) and the max- 
imum Bondi-Hoyle accretion power of ~ 2 x 10 25 W that 
results from the formation of a surrounding HII region 
by the neutron s tar UV radiation, is 10 times too low 
( Blaes et al. 1995 ). The accretion power reached for fast 
(>200-400 km/s), highly magnetized (10 12 G) neutron 
stars with long periods in the intercloud medium (10 -3 
H cm -3 ), though increased by Kelvin- Helmholtz instabil- 
ities in the shocked gas, is also orders of magnitude too 
low (Harding & Leventhal 1992). Isolated accreting black 
holes with masses of 10 M ( |Colpi et al. 1986j ) and 35 
M Q (Dermer 1997) have been proposed, but they would 



be too rare in the Belt: for mass progenitors > 25 M© 



(Timmcs et al. 1996), black holes are 3 to 9 times fewer 
than neutron stars for T indices of —1.1 & —2.0, respec- 
tively. Moreover, the luminosity L^/Ly^ ratios ^> 50 ob- 
served assuming statistically 1 or 2 faint ROSAT source in 
an EGRET error box, are clearly at variance with p-p in- 
teractions in the accretion flow (Colpi et al. 1986), unless 



non-thermal acceleration is advocated (in micro-quasar 
jets?). These luminosity ratios are also at variance with 
standard accreting binary systems. I found no spatial co- 
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incidence with a WR star or the numerous O stars at 
mid latitude despite their highly supersonic winds with 
kinetic powers of 10 28-29 W. Eight pulsars are known 
in 7 rays, 7 bright young ones at kpc distances in the 
Galactic disc, and a faint, older one inside the Gould 
Belt (Geminga). So, the present sample is strongly biased 
to high luminosity and youth. The stability of most of 
the Belt sources ( [Tompkins 1999| , pehrels et al. 2000[ ) is 
consi stent with a pulsar o rigin. Born with large veloci- 
ties ( Lyne fc Lorimer 1994 ), Galactic pulsars rapidly mi- 
grate a way from the plane to mid-lat itudes. The outer gap 
model ( |Yadigaroglu fc Romani 1997 ) for beamed emission 
predicts that 4-5 Galactic pulsars should be detectable at 
|6| > 5° in contrast with the 40 ± 5 sources associated with 
the Belt. Similarly, the wide-beam comptonized polar- 
cap model ( Sturncr fc Dcrmcr 1996 ) predicts 1-2 Galac- 
tic pulsars detectable at |6| > 10° as opposed to 25 ± 
5 sources linked to the Belt. These discrepancies cannot 
be resolved by increasing the Galactic pulsar birth rate 
by more than 30% for fear of overproducing sources at 
low latitude (Yadigaroglu & Romani 1997), nor by using 



larger scale heights or velocities which are not supported 
by the radio data. 



Given the enhanced SN rate in the Belt and its in- 
clined geometry, I propose that the sources associated 
with the Belt be relics of Belt supernovae in the form 
of few Myr old pulsars. Detecting 20 to 40 Belt col- 
lapsed stars as EGRET sources requires the product of 
the beaming fraction AS1/47T and pulsar age be of or- 
der 1-1.5, for instance Att/Air — 0.5 over 2 or 3 Myr. 
AH/Att ~ 0.1 is predic ted for the main polar-cap beam 
( Thompson ct al. 1997 ), and values of 0.1-0.6 are possi- 
blc for the ou ter-gap fan beam depending on pulsar age 
( Romani 199£ ). Yet, one should bear in mind the extreme 
closeness of the Belt objects. The 7-ray luminosity, L 7 , 



scales with the spin-down power, E, as L 7 oc E 1 / 2 over 4 



decades in E ( Thompson et al. 1997 ). Extrapolating from 
the faint Geminga for only half a decade, to E = 10 27 
W => L 7 ~ 6 x 10 25 W over 1 sr, suggests that a pul- 
sar 10 times as old as Geminga, i.e. 3-Myr old, remains 
easily detectable by EGRET out to 500 pc. Furthermore, 
the recorded 7-ray lightcurves show that 10 times fainter 
emission is detected off the main peaks over large phase 
intervals. This side emission from a 3 Myr old pulsar would 
remain detectable by EGRET up to 350 pc, thus largely 
increasing the detection probability An/An. Recent polar 
cap simulations indicate that 4-5 times as many off-beam 
sources as on-beam ones would be detectable above 100 



MeV at a given distance (Harding & Zhang 2001), there- 



fore up to 350 pc for side emission. In this case, the pop- 
ulation of Belt neutron stars born in the last 2-3 Myr 
may account for the Belt 7-ray sources. It may also ex- 
plain the scarcity of bright on-beam-like sources off the 
Galactic plane. The softness of side emission (I7I = 1.8 to 
2.5) may also explain the soft average spectral index \j\ 



= 2.25 ± 0.03 measured for the Belt sources as opposed 
to that of I7I = 1.74 ± 0.02 obtained the 5 on-beam pul- 
sars in the Galactic plane. Preliminary simulations show 
that the Belt spatial signature is preserved over at least 
2 Myr despite rapid migration. So, given our present un- 
derstanding of pulsar 7-ray emission, the hypothesis that 
the Gould Belt sources be powered by pulsars, mostly off- 
beam pulsars, is quite plausible. If true, the Belt pulsars 
would considerably broaden our understanding of photon- 
particle cascades inside their magnetospheres to older and 
lower-luminosity objects at various aspect angles. The ra- 
dio beam being apparently much narrower than the 7-ray 
beam, at least in one dimension, one would expect a ma- 
jority of radio-silent pulsars among the Belt sources. The 
spectral criterion I7I < 2 often adopted to search for 7- 
ray pulsars may not be valid for nearby objects. The next 
generation telescope, GLAST, to be launched in 2005, will 
be able to detect their periodicity in the 7-ray signal, if 
any. As supernova relics, the Belt sources would bring use- 
ful constraints on the initial star mass spectrum at large 
masses and on the abundance of explosive nucleosynthesis 
products in supernova remnants. 

Acknowledgements. With sincere thanks to C. Dermer, N. 
Gehrels, A. Harding, P. Michelson, J. Paul and D. Thompson. 
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Abstract. The nearby, 30 to 40 Myr old, starburst re- 
gion of the Gould Belt has formed numerous massive stars. 
Within its ~300 pc radius, it produces core-collapse su- 
pernovae at an enhanced rate which is shown to be 75 
to 95 Myr -1 kpc -2 , i.e. 3 to 5 times higher than the lo- 
cal Galactic rate, over the past and future few million 
years. A population of persistent, but unidentified, 7-ray 
sources has been recently singled out at medium latitudes 
above 100 MeV. Their distribution across the sky is shown 
to be quite significantly and better correlated with the 
tilted Gould Belt than with other Galactic structures. As 
many as 40 ± 5 sources are statistically associated with 
the Belt at |6| > 5°. It is therefore proposed that these 
sources are part of the Belt and are relics of the Belt su- 
pernovae in the form of million-year old pulsars. Their 
presence stresses how active the local medium is, heated, 
enriched and shaped by multiple recent explosions. 

Key words: Gamma rays: observations - Stars: pulsars: 
general - Stars: supernovae: general 



1. Introduction 

Above 100 MeV, EGRET has discovered ~70 persistent, 
unidentified, sources at medium latitudes that clearly 
differ from the ~40 sources seen close to the Galactic 
plane (Hartman et al. 1999): they are significantly softer, 
fainter, and have a steeper logN-logS function than at 
low latitudes (Gehrels et al. 2000). They have been ten- 
tatively associated with the local interstellar medium 
and the Gould Belt (Grenier 1997; Gehrels et al. 2000; 
Grenier 2000), but their nature remains a mystery. The 
likely candidates present in the Galactic disc, such as pul- 
sars, supernova (SN) remnants, and OB associations, can- 
not account for so many sources off the plane. Interest- 
ingly, the starburst Gould Belt disc is tilted at ~ 20° to 
the Milky Way and close enough for weak sources to be 
detected. Proposing an origin of the sources in the Belt 
however requires one a) to confront their spatial distribu- 
ter^ offprint requests to: I. A. Grenier 



tion with that of the Belt and other likely Galactic struc- 
tures given strong observational biases; b) to evaluate the 
number of likely 7-ray emitters, i.e. compact stars, formed 
in the Belt in the recent past. Both aspects are analyzed 
below before discussing the possible nature of the sources. 

2. Supernova rate in the Gould Belt 

Our Galaxy produces on average 2. 51^4 SNe per cen- 
tury (Tammann et al. 1994), ~85% of which arise from 
the core collapse of a massive star. This value and the 
stellar distribution in the Galaxy imply a local rate of 20 
events Myr -1 kpc -2 , in reasonable agreement with the 
29 progenitors Myr -1 kpc -2 found with masses > 8 M 
within 1 kpc from the Sun (Tammann et al. 1994). This 
region, however, includes the starburst Gould Belt with 
its 300 pc radius. Today, the Belt hosts 432 ± 15 pro- 
genitors with masses > 8 M (Comeron et al. 1994) and 
their maximum lifetime will imply a crude minimum rate 
> 40 collapses Myr -1 kpc -2 in a few tens of Myr. Nu- 
merous explosions have lately occurred from the first gen- 
erations of massive stars formed in the Belt. A recent- 
past rate can be inferred from the current stellar con- 
tent of the Belt as a function of mass, for both short- 
lived and long-lived stars, given a few assumptions: 1) a 
stellar initial mass spectrum dN/dM = dM r_1 with a 
range of indices —2.0 < T < —1.1 as measured in nearby 
OB associations (Scalo 1986; Massey et al. 1995); 2) life- 
time estimates as modelled for stars with solar metallicity 
(Z = 0.02; Schaller et al. 1992; Meynet et al. 1994) and 
interpolated in mass; 3) a constant birth rate for simplic- 
ity; 4) a conservative mass threshold for collapse of 8 M Q . 
Using this formalism, star counts were integrated in var- 
ious spectral bands and scaled to the observed Belt star 
counts (Comeron et al. 1994) to eliminate the unknown 
birth rate and amplitude a. 

The significant increase in birth rate obtained for stars 
<40 Myr confirms the Belt nuclear age and the validity 
of our simple formalism. SN yields of 21.4 ± 0.8, 24.1 ± 
0.9, and 27.2 ± 1.0 per Myr were obtained for T indices 
of —2.0, —1.5, and —1.1, respectively, for a 40 Myr old 
Belt. The quoted error results from the uncertainty in the 
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Fig. 1. visibility map, in Galactic coordinates, from the best 
fit "Iso+Belt" model, the Belt being traced by its young mas- 
sive stars. The colour codes the number of sources detectable 
per bin in the EGRET survey. The 67 persistent, unidentified 
EGRET sources at |6| > 2.5° are marked as white dots which 
closely follow the curved lane of Belt stars. 

observed star counts. These estimates decrease by 15% 
with Belt age from 50 to 30 Myr. Truncating the mass 
spectrum at 60 or 120 M Q , or choosing 10 M for the col- 
lapse threshold, has < 10% impact on the results. Given 
the uncertainties for the Belt age and T index (particu- 
larly at large masses), we infer a frequency of 20 to 27 
SNe per Myr and a rate of 75 to 95 Myr -1 kpc -2 which 
is 3 to 5 times the local Galactic one. This high rate is 
valid for the past few Myr. It is consistent with the power 
of 2.3 x 10 51 erg Myr -1 kpc~ 2 required to maintain the 
local cosmic-ray density (Blandford & Ostriker 1980) for 
a standard SN-to-cosmic-ray energy conversion efficiency 
of a few percent (Drury et al. 1994). It is consistent with 
the presence, within the Belt, of four 0.1-1 Myr radio loops 
(Berkhuijsen 1973) and the Local Bubble. Thirteen radio 
pulsars from the Princeton catalogue are found at high lat- 
itudes with distances <1 kpc and age <2 Myr, but they 
are too few to show a correlation with the inclined Belt 
or with the Galactic plane. The narrow radio beams from 
many more may miss the Earth. 

3. Gamma-ray sources in the Gould Belt 

The distribution on the sky of a subset of the unidenti- 
fied EGRET sources is strongly reminiscent of the tilted 
Gould Belt. The subset includes all persistent sources at 
|6| > 2.5° from the 3rd EGRET catalogue, i.e. those de- 
tected with a significance > 4a (> 5er at \b\ < 10°) using 
the cumulative data from April 1991 to October 1995. 
Sources listed as likely artifacts were discarded. The 67 
sources thus selected are displayed in Fig. 1. 80% of them 
show no time variability, the others are only moderately 
variable (Tompkins 1999). Strong biases due to the in- 
strument exposure and the bright interstellar background 
does influence their apparent spatial distribution. To as- 
sess their correlation with various structures, a maximum- 



likelihood test was applied that takes these biases into 
account (Grenier 1997). An important degree of freedom 
is the unknown luminosity function dN/dL oc L" of the 
parent sources. The flatter the function, the sharper the 
longitude and latitude profiles of any Galactic population 
because more sources remain visible to large distances. 
This effect being quite strong, a was kept as a free param- 
eter over a wide range of values: —2.5 < a < 0. EGRET 
exposure maps above 100 MeV for the 4.5-year survey 
(Hartman et al. 1999), the observed isotropic background 
of (2.1 ± 0.3) x 10~ 5 7 cm -2 s -1 sr -1 , and the interstel- 
lar background model (Hunter et al. 1997) of the EGRET 
survey were used to model the source visibility in 5° by 
5° bins across the sky. 

To allow for possible extragalactic sources among the 
unidentified ones, linear combinations of an isotropic and 
a Galactic component were tested (and nicknamed as in 
Table 1). The choice of Galactic distributions cover the 
various sites likely to display sources at mid latitudes: 
1) sources uniformly distributed in a spherical Galactic 
halo, 20 kpc in radius; 2) sources uniformly distributed in 
longitude in the "local Galactic disc" with any exponen- 
tial or gaussian scale height; 3) sources spread in a "thick 
Galaxy" with a gaussian radial scale length of 9.3 kpc and 
an exponential scale height of 0.4 kpc, typical of radio 
pulsars (Lyne et al. 1985); 4) sources distributed in the 
interstellar medium as mapped by the NH=N(HI+2H 2 ) 
column-densities (Grenier 1997); 5) sources in the Gould 
Belt as traced by the column-densities of stars with spec- 
tral type < B3, fitted from the Yale Bright Star cat- 
alogue data according to the function N*(l,b) oc exp(- 
0.5*sin 2 [6 - &*(;)]/sin 2 <5), with the resulting width S = 
12.0° ±0.4° and median latitude 6*(Z) of the star distribu- 
tion displayed in Fig. 1. Distributions 4) and 5) provide 
independent ways to trace the Belt through its young stars 
as well as its clouds since the latter dominate the inter- 
stellar maps at medium latitudes. 

The reliability of the analysis was checked by applying 
it to the 67 active galactic nuclei discovered by EGRET 
and spread over the sky. For them, an isotropic distribu- 
tion indeed yields the best fit and no Galactic component, 
however large in scale height, is detected: the maximum 
log- likelihood value remains constant within 0.1 for all the 
models tested in Table 1. 67 ± 8 sources are attributed to 
the isotropic component, ± 6 to the local Gal., thick 
Gal., NH, or Belt components. The 9 ± 16 sources at- 
tributed to a halo component are not significant. The lu- 
minosity index of the parent population is a = 1.7 ± 0.3. 

Results obtained for the 67 unidentified sources are dis- 
played in Table 1. Modelled source counts in the isotropic 
and anisotropic components, N- mQ and N alj ^ respectively, 
were summed over the given latitude intervals. The prob- 
ability P a [2|l] is determined from the log-likelihood in- 
crease between two models, with a left free. It is the 
chance probability that a random fluctuation from model 1 
yields as good a fit as model 2. The dramatic improvement 
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Table 1. max-likelihood results and their la errors for selected models against the 67 persistent unidentified sources 
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Fig. 2. latitude profiles of the persistent unidentified EGRET 
sources (crosses), of the best fit "Iso+Belt" model (thick line), 
and of the isotropic contribution to this model (dotted line). 
The "Iso+Local Gal." model (dashed line) yields a significantly 
poorer fit (P a [BELT|LOC] = 3 x 10" 6 ) with systematic devi- 
ations from the data over wide latitude bands. 



in the fit for any Galactic model over a pure isotropic pop- 
ulation shows that over ~50 sources have a Galactic origin. 
The fit improves very significantly from a spherical halo to 
a flatter Galactic distribution, P a [GAL|HAL] = 2 x 1CT 5 , 
and even more so with the tilted geometry of the Belt 
clouds or stars, P a [BELT|HAL] = 2 x 1(T 10 . The good 
quality of the BELT fit is illustrated in Fig. 1 & 2. The 
fact that the LOC and GAL models yield equally good 
fits implies a large fraction of nearby sources: no contrast 
in longitude is detected that would require distances of 
a few kpc to the sources. Observed source counts are in- 
deed equivalent in the centre and anticentre quadrants 
in Fig. 1. Not only are the sources local, the significant 
likelihood increase between the LOC and BELT fits, with 
a chance probability P a [BELT|LOC] = 3 x 1(T 6 , gives 
strong support to their origin in the inclined Belt system. 
On the other hand, the sources at |6| > 5° were shown 
to be distinctly fainter and softer than those at lower lat- 
itudes and a subset of ~20 were pointed out along the 



Belt (Gehrels et al. 2000). Together these findings provide 
compelling evidence that a distinct population of 20 to 40 
EGRET sources belong to the Gould Belt. Based on their 
spatial distribution only, they could be as numerous as 
40 ± 5 at |6| > 5°. The luminosity index of their parent 
population is a = 2.2 ± 0.3. 



4. Nature of the sources 

With luminosities L 7 of 1 to 15 x 10 26 W over 4tt sr 
at 500 pc (for E~ 2 spectra), the sources cannot be un- 
resolved gas clumps irradiated by the local cosmic-ray 
flux: the required mass of ~ 10 4 M Q at 500 pc can- 
not have escaped the radio and IR surveys, even when 
considering radio beam dilution from unresolved clouds 
(Grenier 1997). Nor can they be slow (<20 km/s), old 
neutron stars, accreting gas from a dense cloud: they are 
at least 10 2 ~ 3 times too rare (Grenier 1997) and the max- 
imum Bondi-Hoyle accretion power of ~ 2 x 10 2B W that 
results from the formation of a surrounding HII region 
by the neutron star UV radiation, is 10 times too low 
(Blaes et al. 1995). The accretion power reached for fast 
(>200-400 km/s), highly magnetized (10 12 G) neutron 
stars with long periods in the intercloud medium (10~ 3 
H cm -3 ), though increased by Kelvin-Helmholtz instabil- 
ities in the shocked gas, is also orders of magnitude too 
low (Harding & Leventhal 1992). Isolated accreting black 
holes with masses of 10 M (Colpi et al. 1986) and 35 
M© (Dermer 1997) have been proposed, but they would 
be too rare in the Belt: for mass progenitors > 25 M© 
(Timmes et al. 1996), black holes are 3 to 9 times fewer 
than neutron stars for T indices of —1.1 & —2.0, respec- 
tively. Moreover, the luminosity L 7 /L^ ratios ~3> 50 ob- 
served assuming statistically 1 or 2 faint ROSAT source in 
an EGRET error box, are clearly at variance with p-p in- 
teractions in the accretion flow (Colpi et al. 1986), unless 
non-thermal acceleration is advocated (in micro-quasar 
jets?). These luminosity ratios are also at variance with 
standard accreting binary systems. I found no spatial co- 
incidence with a WR star or the numerous O stars at 
mid latitude despite their highly supersonic winds with 
kinetic powers of 10 28-29 W. Eight pulsars are known 
in 7 rays, 7 bright young ones at kpc distances in the 
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Galactic disc, and a faint, older one inside the Gould 
Belt (Geminga). So, the present sample is strongly biased 
to high luminosity and youth. The stability of most of 
the Belt sources (Tompkins 1999, Gehrels et al. 2000) is 
consistent with a pulsar origin. Born with large veloci- 
ties (Lyne & Lorimer 1994), Galactic pulsars rapidly mi- 
grate away from the plane to mid-latitudes. The outer gap 
model (Yadigaroglu & Romani 1997) for beamed emission 
predicts that 4-5 Galactic pulsars should be detectable at 
|6| > 5° in contrast with the 40 ± 5 sources associated with 
the Belt. Similarly, the wide-beam comptonized polar- 
cap model (Sturner & Dermer 1996) predicts 1-2 Galac- 
tic pulsars detectable at |6| > 10° as opposed to 25 ± 
5 sources linked to the Belt. These discrepancies cannot 
be resolved by increasing the Galactic pulsar birth rate 
by more than 30% for fear of overproducing sources at 
low latitude (Yadigaroglu & Romani 1997), nor by using 
larger scale heights or velocities which are not supported 
by the radio data. 

Given the enhanced SN rate in the Belt and its in- 
clined geometry, I propose that the sources associated 
with the Belt be relics of Belt supernovae in the form 
of few Myr old pulsars. Detecting 20 to 40 Belt col- 
lapsed stars as EGRET sources requires the product of 
the beaming fraction Af2/47r and pulsar age be of or- 
der 1-1.5, for instance Afi/47r = 0.5 over 2 or 3 Myr. 
Af2/47r ~ 0.1 is predicted for the main polar-cap beam 
(Thompson et al. 1997), and values of 0.1-0.6 are possi- 
ble for the outer-gap fan beam depending on pulsar age 
(Romani 1996). Yet, one should bear in mind the extreme 
closeness of the Belt objects. The 7-ray luminosity, L 7 , 
scales with the spin-down power, E, as L 7 oc E 1 /' 2 over 4 
decades in E (Thompson et al. 1997). Extrapolating from 
the faint Geminga for only half a decade, to E = 10 27 
W ==>■ L 7 ~ 6 x 10 2B W over 1 sr, suggests that a pul- 
sar 10 times as old as Geminga, i.e. 3-Myr old, remains 
easily detectable by EGRET out to 500 pc. Furthermore, 
the recorded 7-ray lightcurves show that 10 times fainter 
emission is detected off the main peaks over large phase 
intervals. This side emission from a 3 Myr old pulsar would 
remain detectable by EGRET up to 350 pc, thus largely 
increasing the detection probability Afl/4n. Recent polar 
cap simulations indicate that 4-5 times as many off-beam 
sources as on-beam ones would be detectable above 100 
MeV at a given distance (Harding & Zhang 2001), there- 
fore up to 350 pc for side emission. In this case, the pop- 
ulation of Belt neutron stars born in the last 2-3 Myr 
may account for the Belt 7-ray sources. It may also ex- 
plain the scarcity of bright on-beam-like sources off the 
Galactic plane. The softness of side emission (|7| = 1.8 to 
2.5) may also explain the soft average spectral index [y\ 
= 2.25 ± 0.03 measured for the Belt sources as opposed 
to that of I7I = 1.74 ± 0.02 obtained the 5 on-beam pul- 
sars in the Galactic plane. Preliminary simulations show 
that the Belt spatial signature is preserved over at least 



2 Myr despite rapid migration. So, given our present un- 
derstanding of pulsar 7-ray emission, the hypothesis that 
the Gould Belt sources be powered by pulsars, mostly off- 
beam pulsars, is quite plausible. If true, the Belt pulsars 
would considerably broaden our understanding of photon- 
particle cascades inside their magnetospheres to older and 
lower-luminosity objects at various aspect angles. The ra- 
dio beam being apparently much narrower than the 7-ray 
beam, at least in one dimension, one would expect a ma- 
jority of radio-silent pulsars among the Belt sources. The 
spectral criterion |— ^ | < 2 often adopted to search for 7- 
ray pulsars may not be valid for nearby objects. The next 
generation telescope, GLAST, to be launched in 2005, will 
be able to detect their periodicity in the 7-ray signal, if 
any. As supernova relics, the Belt sources would bring use- 
ful constraints on the initial star mass spectrum at large 
masses and on the abundance of explosive nucleosynthesis 
products in supernova remnants. 
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